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AN  ANALYSIS  OF  HIGH-FREQUENCY  AMBIENT  NOISE 

by 

Lewis  J.  Lloyd  and  Maurice  Daintith 


ABSTRACT 


Because  of  the  high  values  of  absorption  and  bottom-reflection  losses, 
ambient  noise  in  shallow  water  at  frequencies  above  30  kHz  is  considered  to 
result  almost  entirely  from  wind  and  rain  activity  at  the  local  sea 
surface.  The  ambient  noise  is  thus  effectively  produced  by  a  series  of 
small  radiators  distributed  over  a  plane  surface,  with  the  level  and 
directionality  of  the  sources  being  dependent  on  wind  speed,  rainfall,  and 
frequency.  In  order  to  assess  the  form  of  signal  processing  to  use  in  a 
high-frequency  ambient-noise  field,  analytical  expressions  have  been 
developed  for  the  noise  power  in  an  omnidirectional  hydrophone,  for  the 
horizontal  and  vertical  spatial  correlation  functions,  and  for  the  noise 
power  in  an  unshaded  rectangular  array. 


INTRODUCTION 

In  the  early  studies  of  anisotropic  noise  fields  the  effects  of  absorption 
loss  and  bottom  acoustic  interactions  were  neglected  <1,2>.  Omission  of 
the  absorption  loss  implies  that  the  analysis  relates  to  low  frequencies, 
in  which  case  bottom  interactions  should  be  included.  Recent  studies  by 
Kuperman  and  Ingenito  <3>  and  Buckingham  <4>  apply  to  low  frequencies, 
since  bottom  interactions  are  incorporated  but  the  absorption  loss  is 
omitted. 

In  the  study  of  high-frequency  ambient  noise  in  shallow  water  at 
frequencies  above  30  kHz  the  absorption  loss  should  be  included,  since  it 
can  be  an  appreciable  part  of  the  propagation  loss;  for  example,  the 
attenuation  parameter  a  is  equal  to  13  dB/km  at  50  kHz  for  a  temperature  of 
15.5°C.  However,  bottom-reflection  losses  are  of  the  order  of  15  dB  and 
more  per  bounce  for  grazing  angles  between  10°  to  90°,  and  propagation  in 
the  sea-bed  is  negligible,  so  that  bottom  interactions  can  be  neglected. 

In  consequence,  a  good  representation  of  the  high-frequency  ambient-noise 
field  in  shallow  water  is  given  by  assuming  a  number  of  small  independent 
noise  radiators  spread  over  the  sea  surface,  with  there  being  no  coupling 
back  into  the  sea  via  the  sea-bed.  The  level  and  directionality  of  the 
noise  radiators  is  assumed  to  be  a  function  of  the  wind  and  rain  activity 
on  the  local  sea  surface  and  of  frequency.  Because  bottom  interactions  are 
neglected  the  propagation  out  to  ranges  of  about  three  times  the  water 
depth  is  reasonably  well  defined  by  straight-line  ray  paths  and  a 
propagation  loss  given  by  20  log  r  +  or. 
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This  memorandum  extends  the  analysis  of  Cron  and  Sherman  <1>  to  incorporate 

an  absorption  loss  and  provides  expressions  for  the  noise  power  in  an  » 

omnidirectional  hydrophone,  the  noise  level  at  the  output  of  a  conventional 

beamformer,  and  the  spatial  cross-correlation  function.  A  comment  is  also 

made  on  the  likely  variation  of  noise  level  with  frequency.  These  results 

are  of  use  in  assessing  the  performance  of  certain  signal-processing 

techniques  in  high-frequency  shallow-water  ambient  noise. 


1  BASIC  CONSIDERATIONS 

The  model  generally  used  for  low-frequency  ambient  noise  based  on 
theoretical  considerations  <1>  and  experimental  data  <2>  is  that  the  sea 
surface  acts  as  a  series  of  small  independent  noise  sources  radiating 
downwards  into  the  sea.  The  level  and  radiation  pattern  of  the  sources  are 
taken  to  be  a  function  of  wind  speed,  rainfall  and  frequency.  Recent 
studies  have  incorporated  acoustic  interactions  at  the  sea-bed  <3,4>. 

The  small  radiators  are  assumed  to  be  close  to  the  surface  and  to  have  an 

amplitude  radiation  pattern  following  cosmp,  where  p  is  the  angle  measured 
from  the  vertical  and  m  ,  the  surface  noise  directionality  parameter, 
corresponds  to  whether  the  sources  are  monopole,  dipole,  etc.  With  m  equal 
to  zero,  radiation  is  omnidirectional  and  the  sources  are  monopoles; 
with  m  equal  to  1,  they  represent  dipoles  produced  by  specular  reflection 
in  the  sea  surface.  The  dipole  radiation  pattern  of  cosp  is,  however, 
strictly  applicable  only  to  sources  close  to  the  surface  in  terms  of  the 
wavelength,  as  shown  in  Fig.  1. 

At  high  frequencies  (>  30  kHz),  because  acoustic  bottom  interactions  are 
negligible,  the  propagation  out  to  ranges  of  three  times  the  water  depth 
(P  up  to  70°)  is  reasonably  well  defined  by  straight-line  paths.  In 
consequence  the  propagation  loss  is  given  by  20  log  r  +  ar,  where  r  is  the 
slant  range  and  a  the  absorption  coefficient. 

The  use  of  a  cosmp  amplitude  radiation  pattern  at  these  frequencies  may  be 
questionable,  since  we  are  now  dealing  with  small  wavelengths  (<  5  cm)  and 
the  formation  of  dipoles  with  a  separation  of  much  less  than  a  wavelength 
is  not  possible.  However,  it  is  considered  that  some  radiation  pattern  is 
applicable  at  moderate  sea-states  due  to  entrained  air  increasing  the 
near- surface  propagation  loss.  This  radiation  pattern  may  still  be  of  the 

form  cosmp,  since  for  a  given  depth  of  entrained  air  the  path  length 
through  the  layer  increases  as  1/cosp  so  that  the  noise  amplitude  below  the 
layer  may  well  vary  as  cosp,  see  Fig.  2. 

On  the  basis  of  the  above,  the  geometry  of  the  situation  under 
consideration  is  given  in  Fig.  3,  which  shows  a  representative  section  of  a 
small  rectangular  planar  array  located  on  the  sea-bed  with  the  water  depth 
(z)  being  very  much  larger  than  the  hydrophone  spacing  (s).  The  array  is 
shown  tilted  by  an  angle  y  to  the  horizontal  along  the  q  axis  of  the  array. 
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FIG.  1  DIPOLE  RADIATION  PATTERN 
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FIG.  2  RADIATION  THROUGH  A  LAYER  OF  BUBBLES 
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2  NOISE  LEVEL  IN  AN  OMNIDIRECTIONAL  HYDROPHONE 

From  Fig.  3  we  have  that  the  slant  range  r  from  a  small  elemental  area  of 

r  H 

the  sea  surface  to  the  hydrophone  in  the  pth  row  and  the  qth  column  of  the 
array  is  given  as  follows: 


■ ;  I 


rt_  =  (z  -  qs*siny)2  +  R2  , 


where 


„2 

pq 


pq 


R2  =  R2  +  p2s2  +  q2s2*cosy  -  2R(p2s2  +  q2s2*cosy)2  cos(<|)+oO. 
PQ 


Combining  the  above  equations  and  using  (1+x)2  =  l+x/2  if  x  <<1  we  have 


r  =  (z2+R2) 

pq 


1  + 


p2s2+q2s2 


zqs-siny 


2(z2+R2)  (z2+R2) 


Rs 


(z2+R2 ) 


( p - cos(f>  -  q*cosy si n<t>) 


(Eq.  1) 


Thus  we  can  obtain 


rp,q+l  rpq 


(l+2q)s2  zs*siny  Rs*cosysin(|) 


2(z2+R2) 2  (z2+R2) 


(z2+R2) 


Since  s2«z  ,  the  first  term  can  be  dropped,  giving: 


rp  q+1  "  rpq  =  s*cosV*si nP*si n<)>  -  s*sinycosp  . 


(Eq.  la) 


In  a  similar  manner  we  can  obtain 


r  ^  -  r  =  -s*sinB*cos<b. 

p+l,q  pq  h  v 


(Eq.  lb) 


In  the  following  it  is  assumed  that  the  slant  range  (r)  and  the  angle  (p), 
as  measured  to  the  centre  of  the  array,  can  be  used  to  compute  the  noise 
amplitude  at  each  hydrophone  in  the  array.  However,  for  noise  phase  the 
spacing  of  the  hydrophones  is  taken  into  account. 

At  frequency  f  for  the  pq'th  hydrophone,  the  noise  signal  amplitude  in 
a  1  Hz  band  due  to  a  unit  area  on  the  surface  is  given  by: 


sit? 


I  ln 

n*cos  p 


Ar 


fpq 


-j(wt-kr) 


(Eq.  2) 


5  -i-: 

'  4* 

,  ^  ■■ 
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where  n*cosm3  -  noise  signal  radiated  per  unit  area  of  the  surface 
in  a  1  Hz  band, 

A  =  absorption  coefficient  for  signal  amplitude  expressed 

in  exponential  form, 

k  =  wavenumber  =  . 

A 


Then  the  noise  power  from  the  total  surface  is  given  by: 


“>  2  n 

S  J 


o  o 


n 

fpq 


nfpq  R  ^  dR 


(Eq.  2a) 


where  E{  }  denotes  the  expected  value  and  the  asterisk  (*)  indicates  the 
complex  conjugate.  This  gives,  via  Eq.  2  and  Fig.  3, 


n  27t 


n/ 2  „  . 

J  cos  mp*tanp*e  z 


seep 


dp, 


(Eq.  2b) 


where  A  -  a/4.343,  a  being  the  absorption  coefficient  in  dB/m. 


Equation  2b  can  be  readily  integrated  by  using  the  substitution  cosp  =  1/t, 
giving: 

"fpq  =  E2m+1(AZ)  '  (Eq-  2c) 

where  E^  stands  for  the  exponential  integral  of  order  k. 

If  the  absorption  loss  is  omitted,  the  exponential  term  in  Eq.  2b  would  not 
be  present  and: 


(no  absorption  loss)  = 


n2n 

2m 


(Eq.  2d) 


This  suggests  that  at  low  frequencies,  with  no  bottom  interactions,  the 
noise  level  in  an  omnidirectional  hydrophone  (Eq.  2d)  is  independent  of  the 
depth  of  the  hydrophone.  In  our  high-frequency  case  (Eq.  2c),  the  noise 

power  is  a  function  of depth  and  of  absorption  coefficient.  Thus  the 

high-frequency  value  of  N|  has  two  components  that  vary  with  frequency: 

the  value  of  n2,  which  for  a  given  sea  condition  decreases  with  increasing 
frequency,  and  the  exponential  integral  E2m+1(Az),  which  decreases  with 

increasing  frequency  due  to  the  increase  in  A. 

Figure  4  plots  values  of  in  decibels  with  respect  to  the  level  at 

10  kHz  for  three  typical  shallow-water  depths  of  150,  300  and  450  m  for 
frequencies  from  10  to  100  kHz.  The  graphs  are  for  values  of  m  =  0  and 
m  =  3,  the  extreme  values.  The  value  of  m  =  0  is  likely  to  apply  at  low 
frequencies  and  for  only  small  activity  on  the  sea  surface,  that  of  m  =  3 
for  high  surface  activity  and  high  frequencies.  Superimposed  on  these 
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©2  =  150  m 


m  =  0 


A  z  =  300  m 
□  z  =  450  m 


FIG .  4  VARIATION  OF  E 0  ,  MZJ  WITH  FREQUENCY 


1 


graphs  is  a  line  showing  the  generally  accepted  slope  for  high-frequency 
ambient  noise  obtained  by  linear  extrapolation  from  lower  frequencies.  It 
is  seen  that  the  E2m+^(Az)  curves  show  a  greater  fall -off  above  40  kHz  than 

given  by  the  linear  extrapolation,  particularly  for  low  sea-state  (m  =  0) 
and  the  greater  water  depth  (>  300  m).  Because  Eqs.  2c  and  2d  represent 
spectrum  levels,  the  total  noise  power  in  a  band  form  fx  to  f2  is  obtained 
by  integrating  over  the  band: 

_  f  2 

NWpq  =  {  *  "  WAz)df  • 

T  1 

where  band  W  =  f2-fl. 


3  SPATIAL  CORRELATION  FUNCTION 

For  a  single  frequency  (f)  in  our  band  of  interest,  the  normalized  spatial 
correlation  function  is  given  by: 

E{n(f,x1,y,z)*n)k(f  ,x2,y,z)} 

p(x,f)  =  - r-  ,  (Eq.  3) 

{E[n2(f  ,x1,y,z)]E{n2(f  ,x2,y,z)]}  2, 


where  E{  }  refers  to  the  expected  value  of  the  expression,  n(f,x,y,z) 
represents  the  noise  signal  at  frequency  f  and  position  (x,y,z),  and  it  is 
assumed  that  the  noise  signals  have  zero  mean  values. 

In  terms  of  our  geometry  and  model,  and  following  the  form  of  Eqs.  2  and 
2a,  Eq.  3  can  be  written  as: 

1  n/2  2n 

P(s,f)  =  —  f  f  E  nfpq  -n*p  q+1  R  d*  dp,  (Eq.  3a) 

fpq  '  ’ 

where  has  been  taken  equal  to  N?  on  the  basis  of  our  earlier 

assumptions. 


Thus  the  spatial  correlation  function  for  the  noise  from  the  total  surface, 
following  the  form  of  Eqs.  2b  and  2c,  is  given  by: 


p(s,f) 


(n/2  2 n  n2  cos^mp  -Ar 

E  J  S  -  e  • 

'  o  o  2r2 


j  k( r  ^T-r  ) 

,  P,q+1  pq 


R  d(j)  dp 


(Eq.  3b) 
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which  leads  to 


p(s,f)  =  - J  cos2nrlp*sinP'e~Az*Sec^  f  cos(ks*cosY*sinp*sin<|> 


2Nlnn 

fpq 


ks'SinY’Cosp)  d(J>  dp. 


In  the  integral 
we  have: 

p(s,f)  .=  - 


over  <t>  the  integral 


n/2 
f  cos 


2m- 1 


of  sin(ks*cosY*sinp*sin<t>)  is  zero;  thus 
p*sinp  e  Az  secP  cos(ks’COSp*sinY) 


E2m+l(flz)  0 


•J  (ks'cosvsinp)  dp. 

0  (Eg.  3c) 


A  solution  to  this  integral  is  given  in  Appendix  A.  For  a  horizontal 
separation  (y  =  0),  the  following  expression  obtains: 


p(x,f) 


00  (^ks)^ 

I  -  J  (ks)E 

H 


2m+l+2p 


(Az) 


E2m+l(Az)  M=0 
where  k  =  wavenumber  and  s  =  hydrophone  spacing. 


(Eq.  3d) 


For  a  vertical  separation  (y  =  90°)  we  have: 

1 _  “  (-k^)1 


p(z,f)  = 


E2m+l(Az)  1=0  (2i)! 


E2m+l+2i<‘Az)‘ 


(Eq.  3e) 


The  corresponding  equations  by  Cron  and  Sherman  <1>  for  the  case  when 
absorption  loss  is  neglected  are 


P'(x) 


2mm!  Jm(ks) 
(ks)m 


(Eq.  3f ) 


p 1 (z)  =  2m  J  X2m_1  cos(ksX)  dX. 


(Eq.  3g) 


As  can  be  seen  from  Appendix  A,  it  is  not  easy  to  show  directly  that  when  A 
is  set  equal  to  zero  Eqs.  3d  and  3e  are  the  same  as  Eqs.  3f  and  3g 
respectively.  However,  it  can  be  shown  indirectly  by  comparing  the 
numerical  values  obtained  from  the  two  pairs  of  equations. 

This  comparison  is  done  for  different  values  of  m  in  Fig.  5  for  p(x)  and  in 
Fig.  6  for  p(z),  where  the  solid  line  represents  the  values  from  the 
present  analysis  with  a  very  small  value  of  A  (A  is  set  equal  to  0.0001, 
which  represents  0.43  dB/km),  and  the  dashed  lines  represent  the  values 
from  Cron  and  Sherman  <1>. 


9 


SACLANTCEN  SM-159 


HORIZONTAL  CROSS-CORRELATION  FUNCTION 
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FIG .  5  COMPARISON  OF  p (x)  WITH  LOW  FREQUENCY  ANALYSIS 


VERTICAL  CROSS-CORRELATION  FUNCTION 
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FIG.  6  COMPARISON  OF  p (z)  WITH  LOW  FREQUENCY  ANALYSIS 
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It  is  seen  that  good  agreement  exists  with  p'(x)  for  all  values  of  m,  and 
also  with  p'(z)  except  for  the  case  of  m  equal  to  zero.  In  this  instance, 
it  is  seen  that  the  Cron  &  Sherman  expression  gives  a  value  of  1  for  p(z), 
independent  of  the  normalized  element  spacing  s/X,  whereas  when  an 
absorption  loss  is  included  the  present  analysis  gives  a  reducing  value  of 
p(z)  as  s/X  is  increased.  As  will  be  seen  in  the  next  figures,  the  value 
of  p(z)  becomes  oscillatory  about  zero  for  m  equal  to  zero  when  the 
absorption  loss  is  greater  than  a  value  corresponding  to  20  kHz  and  for 
water  depths  greater  than  150  m. 

Value  of  p(x)  and  p(z)  are  plotted  in  Figs.  7a,  7b,  7c  and  7d  and  in 
Figs.  8a,  8h,  8c  and  8d  respectively.  The  graphs  are  for  the  three  sample 
depths  of  150,  300  and  450  m,  frequencies  of  20,  50  and  100  kHz, 
and  m  equal  to  0,  1,  2  and  3. 

The  effects  of  the  three  variables  of  m,  frequency,  and  water  depth  shown 
by  Figs.  7  and  8  are  plausible  and  consistent  with  the  previous  deep-water 
low-frequency  analyses  <1,2,3>.  The  values  of  s/X  at  which  p(x)  and  p(z) 
have  their  first  zeros  are,  for  m  equal  to  1,  a  water  depth  of  450  m,  and  a 
frequency  of  20  kHz,  0.6  and  0.3  respectively;  the  values  at  low  frequency 
<1,2,3>  are  0.6  and  0.35,  see  Figs.  5  and  6.  The  effect  on  p(z)  and  p(z) 
of  increasing  m,  that  is  of  making  the  field  more  directional,  is  also 
consistent  with  the  low-frequency  values.  Increasing  the  water  depth  will 
increase  the  propagation  loss  from  the  surface  along  a  given  angle  of 
incidence  at  the  sea  bed,  which  will  make  the  field  more  directional  and 
thus  have  a  similar  effect  to  increasing  m.  A  similar  argument  applies  to 
an  increase  in  frequency  (increases  in  absorption  coefficient).  In 
summary,  therefore,  the  curve  for  p(x)  is  effectively  stretched  out  by  an 
increase  in  m,  water  depth,  or  frequency,  whereas  the  opposite  applies  to 

p(z). 


4  NOISE  LEVEL  IN  A  BEAM 

The  analysis  given  above  can  be  extended  to  give  the  noise  power  in  a  beam 
produced  by  a  conventional  beamformer,  that  is,  a  beam  formed  by  summing 
the  phase-shifted,  but  unweighted,  hydrophone  outputs. 

The  array,  with  dimensions  small  compared  with  the  water  depth,  is  assumed 
to  be  tilted  at  an  angle  y  to  the  horizontal  in  the  q  direction,  and  the 
steering  is  taken  to  be  along  the  q  axis  only,  so  that  the  beam  is  steered 
in  the  same  plane  as  the  array  is  tilted.  The  beamforming  can  thus  be 
represented  either  by  summing  along  the  q  axis  after  the  appropriate  delays 
and  then  summing  across  p,  or  by  summing  across  p  without  steer  delays  and 
then  summing  along  q  with  the  appropriate  delays.  The  first  sequence  will 
be  used  in  the  following. 

Let  the  steer  angle  measured  from  the  vertical  be  6,  then  the  phase  shifts 
required  along  the  q  axis  are  qk*sin0,  where  q  varies  from  -(Q-l)/2  to 
+(Q-l)/2  or  from  0  to  Q-l. 

For  the  wavefront  of  the  noise  from  a  small  elemental  area  of  the  sea 
surface  a  reference  time  is  taken  at  one  end  of  the  array  (the  'zero1  end 
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FIG.  7a  HORIZONTAL  CROSS-CORRELATION  FUNCTION  (M  =  O) 
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HORIZONTAL  CROSS-CORRELATION  FUNCTION 


p(X) 


E\V 

E 

E 

r  2 

E 

3 

E 

E 

i 

ILLL 

1LLL 

lLLLL 

JILL 

LULL 

JLLL I 

JILL 

JJ-LL 

ILLL 

ILLL 

iiLL 

JUl 

WATER  DEPTH 
1  =  150m 
2=  300m 
3=  450m 


20  kHz 


.00  ,2S  .BO  .75  1.00  1.26  l.SO  t .76  2.00  2.2S  2.50  275  3.00 

/•I 


s/x 


=> 

\ 

2 

> 

3 

\ 

E 

E 

E 

ILLL 

ILLL 

JLLJ1 

JILL 

II Li 

JILL 

IILL 

JILL 

IILL 

IILL 

Jill 

ILLL 

50  kHz 


.00  .25  .50  .75  t  - 00  1  .25  1.50  1.75  2.00  2.25  2. 50  2-75  3  00 


E 

r  2 

E 

3 

E 

E 

E 

E 

Till 

mi 

mi 

JUI 

JUI 

ILLL 

IILL 

IILL 

JUL 

INI 

mi 

ILLL 

100  kHz 


.00  .25  .50  .75  1.00  t.26  l.SO  1.75  2.00  2.26  2.50  2.75  3.00 


s/X 


FIG.  7b  HORIZONTAL  CROSS-CORRELATION  FUNCTION  (M  =  1) 


13 


SACLANTCEN  SM-159 


HORIZONTAL  CROSS-CORRELATION  FUNCTION 
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HORIZONTAL  CROSS-CORRELATION  FUNCTION 
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VERTICAL  CROSS-CORRELATION  FUNCTION 
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FIG.  8a  VERTICAL  CROSS-CORRELATION  FUNCTION  (M  =  O) 
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VERTICAL  CROSS-CORRELATION  FUNCTION 
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FIG .  8b  VERTICAL  CROSS -CORRELATION  FUNCTION  (M  =  1) 
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VERTICAL  CROSS-CORRELATION  FUNCTION 
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FIG.  8c  VERTICAL  CROSS -CORRELATION  FUNCTION  (M  =  2) 
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VERTICAL  CROSS-CORRELATION  FUNCTION 
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FIG.  8d  VERTICAL  CROSS -CORRELATION  FUNCTION  (M  =  3) 
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in  the  above  summation).  Thus  the  incremental  delays  of  the  wavefront 

along  the  q  axis  can  be  obtained  from  the  increment  in  the  slant  range 

(r  )  along  the  q  axis,  i.e.  by  r  .. ,  -  r  . 

pq  p,q+i  pq 

The  noise  amplitude  along  the  q  axis  in  the  pth  column  due  to  a  small 
elemental  area  of  the  surface  is  then: 


„  n'C0S  15  -ju>t  ‘J;1  -MK'51'1'8  -i^rp,qn'rpq)  'flR 

N  (0)  =  -  e  J  I  e  e  e 

9  r  q=o 


with  the  value  of  0,  and  that  of  r  associated  with  noise  amplitude,  taken 
as  independent  of  p  or  q  since  the  array  dimensions  are  small  compared  with 
the  water  depth. 


The  sum  of  N  (0)  across  the  p  axis  can  similarly  be  obtained  by  taking 

q 

increments  of  r  ^  -  r  .  From  Eqs.  la  and  lb,  we  have 

p+l,q  pq 

r  +,  -  r  =  s(cosy'sinp*sin<t>  -  sinycosp) 
p ,  q  x  pq 

and 


rp+i>q 


r  =  -s*sin|3*cos<t>. 

pq 


Then  the  total  noise  amplitude  across  the  full  array  is  given  by: 


(0) 


PI1  e-jpks*sinp-cos<t>  N  m 

p=0  « 


(Eq.  4) 


That  is: 

N  (0)  =  n  cosmp  £  ‘ju,t_Ar  e-jPks*sinP*cos<t> 

Pq  r  p=0 

e-jks(sinp-cosysin())-sinycosp)  e~jks*sin0 

q=0 


=  n  cosmp*-  ■ju,t‘Ar  Pi1  e‘jPx  9i1 
r  p=0  q=0 


where 

x  =  ks*sinp*cos<}) 

and 

y  =  ks(sinp*cosysin(()  -  sinycosp  +  sin0). 


(Eq.  4a) 

(Eq.  4b) 
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Si  nee 

P-1 
2  e 

p=0 


“JPX  _ 


=  e 


•X/r.  -v  •  Px 
-j^CP-l)  sin 


sin 


Q-l  . 

I  e 
q~o 


jqy  _ 


-j|(Q-l)  sin  ^ 


sin 


(Eq.  4c) 


(Eq.  4d) 


Thus  we  have': 


.x,n  •  Px 
•  +.  a  -Jo(P-l)  sin  7T- 
tl  m„  e  .iwt-Ar  J2  2 

N  (0)  =n*cos  B  -  e  - 

pq  r  x 

sin  | 


e-j|(Q-l)  sin  f- 


si  n 


The  total  noise  power  in  the  beam  is  then  given  by: 
(»  2  n 

N^(0)  =  E  /  J  N  (0) »N*  (0)R  d<t>  dR 
PQ  PQ 


o  o 


n/Z  2 n 


=  n2  P2Q2  f  f  cos^m  1B  sinB‘6  ^ Z  SeC ^  [f(P)*F(Q)]2  d<()  dB  , 

(Eq.  4e) 


o  o 


where  fm  =  sintePks-sinp-cosfr) 

'  PsinCJgks'sinB’Coscj)) 

and  F(0)  =  s^n[3sQ*<s(s^nP,cosy‘s''n<i)  ~  sinycosB  +  sinB)] 
Q’sint^ksCsinB'Cosysin^  -  sinycosB  +  sin0)] 


Equation  4e  is  evaluated  in  Appendix  B  to  give 

oo  oo  p-i  Q-i 

W(n,  _  2nn2  2  2  2  2  A  B.  (-)i  (\/2)n  Jn(V) 

N  (6)  -  pSp?  i=0  n=0  r=0  t=0  r  1  ”HT 

C°=P'?21  E2m+2n+2i+l(Az)  +  E2m+2n+2i+2^z)  ’ 


where 


A  =  P 
o 


B0  =  Q 


Ar  =  2(P-r)  (r  >  0) 
X  =  ks  7  r2+t2cos2y 


(2i+l)! 


Bt  =  2(Q- t)  (t  >  0) 
p  =  kst-sin0  q  =  kst*siny 
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This  is  a  fairly  complex  expression  and,  of  course,  applies  only  to  an 
unshaded  array.  A  closed-form  solution  for  a  shaded  array  would  be 
achievable  if  the  shading  functions  allow  summations  to  be  made,  as  in 
Eqs.  4c  and  4d.  If  the  shading  function  is  known,  then  the  total  power  in 
a  beam  is  probably  best  evaluated  using  numerical  integration. 

For  a  horizontal  array,  in  which  y  =  0,  Eq.  4f  becomes 


N*(6)h  = 


Znn2 


P2Q2  n=0  r=0  t=0 


»  P-1  Q-l  (Ai/2) 

z  i  ^  J,w 

(Eq.  4g) 


n! 


where  Ax  =  ksV  r2+t2  , 
p  =  kst*sin0  . 


In  the  case  of  a  vertical  array,  in  which  y  =  n/ 2,  the  quadruple  summation 
remains  as  in  Eq.  4f,  but  now 

A2  =  ksr  q  =  kst 

(A2/2)n 

and  there  are  only  P  terms  of  the  form  -  Jn(A2). 


CONCLUSIONS  AND  RECOMMENDATIONS 


If  the  model  of  high-frequency  ambient  noise  in  shallow  water  proposed  and 
examined  in  this  study  is  reasonable,  then: 

a.  Care  should  be  taken  in  using  noise  levels  obtained  by  linear 
extrapolation  from  the  lower  frequencies.  The  study  suggests  that  the 
noise  power  in  an  omnidirectional  hydrophone  (Eq.  2c  and  Fig.  3)  falls  off 
more  rapidly  with  frequency  above  40  kHz,  particularly  for  low  sea  states 
(m  equal  to  zero)  and  for  greater  water  depth  (z  =  150  m),  than  is 
predicted  by  linear  extrapolation  from  the  lower  frequencies. 

b.  When  designing  an  array  for  use  near  the  sea  bed  at  frequencies 
above  20  kHz,  the  value  of  hydrophone  separation  for  the  first  zero  in  p(x) 
or  p(z)  may  be  a  function  of  sea  state,  frequency,  and  water  depth.  From 
the  results  in  Figs.  7  and  8  it  would  appear  necessary  to  space  the 
elements  by  more  than  A/2  for  horizontal  arrays  and  less  than  A/2  for 
vertical  arrays,  depending  on  the  value  of  m  that  is  applicable,  the 
frequency  in  use,  and  the  water  depth.  For  example,  for  m  equal  to  1,  a 
frequency  of  50  kHz,  and  a  water  depth  of  300  m,  the  optimum  spacing  for  a 
horizontal  array  is  0.75A  and  a  vertical  array  0.3A. 

c.  The  expression  for  the  total  ambient  noise  power  in  an  unshaded 
array  is  complex  (see  Eqs.  4f,  4g  and  4h).  For  a  shaded  array,  it  is 
possible  to  obtain  a  closed-form  solution  only  for  those  shading  functions 
that  allow  the  summations  in  Eqs.  4c  and  4d.  For  other  shading  functions, 
however,  a  solution  could  be  obtained  using  numerical  integration. 
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Experimental  data  are  required  to  verify  the  model  and  the  resulting 
analysis.  Measurements  are  proposed  as  follows: 

1)  Noise  Levels 

The  noise  level  in  an  omnidirectional  hydrophone  placed  on  or  near  the 
sea-bed  should  be  measured  as  a  function  of  sea-state  (wind  speed), 
frequency,  and  water  depth  for  comparison  with  Eq.  2c  and  Fig.  3  and  to 
give  estimates  of  the  value  of  n2. 


2)  Nojse  Directionality 

Measurements  should  be  made  of  the  high-frequency  noise  level  in  shallow 
water  using  a  high-resolution  array  (narrow  mainlobe  and  low  sidelobe 
levels).  The  array  should  be  situated  on  or  near  the  sea-bed  and  capable 
of  being  tilted  to  various  angles  away  from  the  vertical.  Preferably,  for 
a  given  sea-state,  simultaneous  measurements  should  be  made  at  two 
different  depression  angles  e.g.  one  array  at  zero  depression  (looking 
vertical),  the  other  at  30°  depression.  A  sampling  technique  could  be 
developed  to  measure  the  noise  levels  for  a  given  sea-state  over  the  whole 
vertical  arc;  for  example,  with  two  arrays,  sample  vertically  and  at  30° 
depression  for  5  minutes,  then  at  10°  and  40°  for  5  minutes,  then  at  20° 
and  50°,  etc. 


3)  Spatial  Correlation 

It  is  recommended  that  consideration  be  given  to  the  construction  of  a 
device  to  measure  p(x)  and  p(z)  over  the  frequency  range  of  30  to  100  kHz. 
This  would  require  the  use  of  very  small  ball  hydrophones,  say  0.2  cm 
diameter,  spaced  on  horizontal  and  vertical  axes  in  such  a  manner  as  to 
provide  adequate  samples  of  p(x)  and  p(z)  out  to  values  of  s/X  =  2.  This 
could  be  done  with  24  hydrophones  (12  on  each  axis)  with  a  separation  of 
0.5  cm  (0.167X). 


4)  Noise  Power  in  a  Beam 

It  is  recommended  that  the  outputs  of  the  small  hydrophones  outlined  above 
be  summed  to  form  simple,  unsteered  beams,  with  or  without  shading,  to 
provide  experimental  results  for  comparison  with  Eqs.  4g  and  4h.  The 
elements  could  be  tilted  to  evaluate  the  effect  of  y,  the  tilt  angle. 
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APPENDIX  A 

EVALUATION  OF  THE  BESSEL  TRIGONOMETRIC  INTEGRAL  FOR  THE 

SPATIAL  CROSS-CORRELATION  FUNCTION 


In  Eq.  3c  of  the  main  text,  the  spatial  cross-correlation  function  includes 
the  integral 

n /2  2m-1 

Ii  =  J  cos  p*sinp*exp(-Az  secp)cos(ks*cosp*siny)Jo(ks*cosysinp)  dp, 

0  (Eq.  A. 1) 

where  m  is  a  positive  integer,  and  Az  is  a  positive  constant. 


As  this  integral  does  not  seem  to  have  been  evaluated  before  (since  it  does 
not  appear  in  collections  of  integrals  such  as  <A.1>,  the  development  of  a 
form  suitable  for  numerical  calculations  is  presented  here.  For  the  sake 
of  completeness,  the  result  is  given  for  y  arbitrary,  although  from  a 
practical  point  of  view  most  interest  centres  on  the  limiting  cases 
of  y  =  0  and  y  =  n/2. 

An  interesting  outcome  of  this  analysis  is  the  discovery  of  two  identities, 
which  may  be  of  use  elsewhere. 


A.  1  Evaluation  of  the  Integral 

Expanding  cos(ks*cosp*siny)  and  J  (ks-cosysinp)  as  infinite  power  series 
in  cosp  and  sinp  respectively, 

00  oo  n/2  1  (-k2s2sin2ycos2p)^ 

Ii  =  I  I  J  cos  p*sinp*exp(-Az*secp)  - - 

i=0  j=0  o 

(~k2s2cos2y  sin2P/4)J  dp  . 
j!  j! 


The  term  (sin2p)^  has  a  binomial  expansion  in  cos2p,  i.e.  (l-cos2p)^. 
Doing  this,  and  collecting  like  terms  together 


Ii 


00  oo  j  n/2 

II2S 

i=0  j=0  t=0  o 


(-) 


i+2j-t(k2s2) 


1+Jsin2lY(cosy/2)2j» 

(2i )!  j! t! ( j-t)! 


cosp2"1*21*2^1 


sinp  exp(-Az-secp)  dp 


(Eq.  A. 2) 
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ti/2 


N 

But  in  J  cos  p*sinp*exp(-Az*secp)  dp,  putting  sec 
o 

becomes 


=  t  this  integral 


r  exp(-Azt)  .. 

{  tN+2  az 

=  the  exponential  integral  EN+2(Az).  (Eq.  A. 3) 


Applying  this  to  Eq.  A. 2,  gives 

li  =  “  "  1  (-)1't(k2s^)itJs)n2tv(cosv/2)2J 
i=0  j=0  t=0  C2i)!  j!t!  (j-t)! 

Write  t  =  j-n  and  sum  over  j  and  n.  Since  t  >  0  and  j  >  n, 


Ji=  *  “  "  (-)Un'J(t=2s2)itJ  s1n21v(cosv/2)2j  E2„2H2n+1 

i=0  n=0  j=n 


(Az) 


(2i )!  j!  (j-n)!  n! 


(The  change  in  order  of  summation  is  justified  since  it  is  easily  seen  that 
the  triple  summation  is  absolutely  convergent.) 


Now  put  j  =  n  +  r 

00  00  03  (“)1~r(k2s2)1+n+r  sin2ly(cosy/2)2n+2r  E 


I  i=I  I  I 
i=0  n=0  r=0 


2m+2i+2n+l 


(Az) 


(2i )!  n!  r!  (n+r)! 


But 


Jn(z)  =  (z/2)  I 


»  (-z2/4)’ 


r=0  r!  (n+r)! 

and  so  the  expression  becomes,  on  summing  over  r, 
00  00  (-k2s2ysin2y)1  (ks*cosy/2)n 


Ii  =  2  2 

i=0  n=0 


(2i )! 


n! 


■Vks-cosY)  E2m.21+2n+1(Az) 

(Eq.  A. 4) 


Equation  A. 4  is  the  required  closed  form.  As  written,  it  is  a 
Fourier-Bessel  expansion 


xi  =  ^  A  J  (ks-cosy) 
n=0  n  n 
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in  which 


(ks*cosy/2)n  oo  (-k2s2*sin2y)1 

A-  = - 7? -  " 

n’  i=0 


(2i )! 


E2m+2i+2n+l(Az)  ' 


For  computing,  however,  it  may  be  better  to  expand  in  a  series  of 
exponential  integrals.  Writing  i  +  n  =  p  in  Eq.  A. 4,  Ix  becomes 


*  Bp  E2m+2p+l(Az)’ 


where 


B 


P 


P  (-k2s2*sin2y)*3  n  (ks*cosy/2)n 

„Z=0  S 


Jn(ks-cosy) 


The  coefficient  B 
advantage.  ^ 


contains  only  a  finite  number  of  terms,  which  may  be  an 


This  general  form  will  not  be  examined  in  further  detail;  as  stated  in  the 
main  text,  the  significant  cases  are  those  for  which  y  =  0  and  y  =  n/2. 
Each  of  these  leads  to  a  significant  simplification  in  the  expression  for 
Ii- 


A. 2  Significant  Cases 
A. 2.1  Horizontal  Array  (y  =0) 

In  this  case,  from  Eq.  A. 4,  in  the  summation  over  i  the  only  non-zero  term 
is  that  for  which  i  =  0,  giving 

I,(v  =  0)  =  l  Jn(ks)  E2m+2n+1(Az).  (Eq.  A.5) 

n=0 


This  is  evidently  a  convergent  series,  since  for  n  >  0  the  moduli  of  both 
the  Bessel  function  and  the  exponential  integral  are  less  than  unity,  so 

(ks/2)n  i.e.  less 
n! 


that  the  moduli  of  successive  terms  in  n  are  less  than 


than  the  corresponding  terms  of  the  absolutely  convergent  series  for 
exp(ks/2). 


Numerical  Computation 


The  following  remarks  will  be  superfluous  to  those  who  are  practiced  in 
numerical  techniques,  but  may  be  useful  to  those  who  are  not  specialists, 
and  who  (as  has  happened  so  often  to  the  author)  can  easily  produce  invalid 
answers. 

In  summing  a  convergent  series  of  the  form  uo+u1+y2+- • • +un- . .  it  is 

customary  to  select  a  desired  accuracy  e  (e.g.  e  =  10  for  accuracy  to  the 
third  decimal  place),  and  to  terminate  the  series  when  |un|<  e •  This 
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assumes,  however  (among  other  things)  that  |un|  decreases  monotonically 

with  n.  This  condition  is  not  necessarily  true  for  the  series  of  Eq.  A. 5. 

First,  Jn(ks)  is  an  oscillatory  term;  i.e.  ,  it  passes  through  zero. 

If,  for  example,  ks  happened  to  be  equal  to,  or  close  to,  a  value 

corresponding  to  a  zero  of  Jn(ks),  the  naive  application  of  the 

rule  | un| c  s  would  stop  the  series  at  n  =  1,  irrespective  of  the  error. 

Secondly,  the  criterion  |un|<  6  is,  in  any  case,  only  fully  justified  if 

the  series  consists  of  terms  alternatively  positive  and  negative.  (For 

00  1 

example,  the  series  I,|  —  diverges  to  infinity,  and  stopping  at  n  =  100 
does  not  give  an  answer  correct  to  two  decimal  places.) 

These  difficulties  may  be  overcome  by  finding  an  upper  bound  to  the  residue 
obtained  by  terminating  the  summation  at  the  n'th  term. 

Since,  for  the  purpose  in  the  main  text,  the  required  quantity  is 
normalized  by  dividing  by  the  value  of  I,  for  ks  =  0,  viz  E0  ,(Az),  we 
have  2m+1 


where 


u 


n 


Hence  the  modulus  of  the  residue  R  on  stopping  at  the  nth  term  is 


00 


R  =  I  u  < 
_ n  - 


r=,,  ■  x 


Now 


Jn(ks)  <  1  for  n  >  0 


and 


^2m+2n+l^z)  <  1  for  n  >  0. 


Hence 


(ks/2)n 


n! 
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and 


I  ( ks/2)n+1  +  (ks/2)n+2  + 

K|  <  (n+1)!  (n+2)! 


(ks/2) 


n+1 


(n+1)! 


(ks/2)  (ks/2)2 


1  + 


n+2 


(n+3) (n+2) 


(ks/2) 


n+1 


(n+1)! 


(ks/2)  (ks/2)2 


1  + 


n+2 


(n+2)5 


(ks/2) 


n+1 


(n+1)! 


(ks/2) 

(n+2) 


A  sufficient  condition  for  terminating  the  summation  is  thus 
.n+1 


(ks/2)' 

(n+1)! 


1  - 


(ks/2) 

n+2 


<  £  . 


(Eq.  / 


This  does  not  give  unmanageable  quantities.  For  example,  if  ks  =  6,  n 

-3  _  5 

gives  an  error  <  0.3  x  10  and  n  =  13,  an  error  <  7  x  10  . 


Two  identities 


If  Az  =  0  in  Eq.  A. 5 


Ia(V  =  0,  Az  =  0)  = 


(since  E  (0)  =  — =•  . 

P  p-1  ). 


00  (ks/2)  Jn(ks) 

I 

n=0 


n!  2(m+n) 


Putting  -y  =  Az  =  0  in  the  original  equation  (Eq.  A.l)  gives 

tt/2  „  , 

Ii(Y  =0,  Az  -  0)  =  2  f  co$*m  p*sinp*Jo(ks*sinp)  dp. 

o 


•  6) 


=  12 
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But  this  is  a  known  integral  (see  <A2>  for  example),  which  integrates  to 
Ii(Y  =  0,  Az  =  0)  =  2m"1(m-l)!  Jm(ks)(ks)"m 

and  so  we  have  the  identity 

(ks/2)m  «  (ks/2)nJ  (ks) 

Jra(kS)  S  TiFijr  Z-  (m-n)-n!  (Eq.  A.7) 

n=0  J 

which,  oddly  enough,  expresses  Jm(ks)  as  a  Fourier-Bessel  series. 


From  Eq.  A.7  an  algebraic  identity  can  be  deduced.  Expanding  the  Bessel 
functions  as  power  series  in  ks 


(ks/2) 


m 


oo  (-k2s2/4;  (ks/2) 


m 


oo  (ks/2)2n  (-k2s2/4)^ 


I  I 


r=0  r!(m+r)!  (m-1)!  n=0  t=0  (m+n)‘n!  t!  (n+t)! 


or 


oo  (-)r(k2s2/4)r 


oo  oo  (-)t(k2s2/4)n+t 
1  I 


r=0  r!(m+r)!  (m-1)!  n^0  t=0  (m+n)  •  n!  t!  (n+t)! 


writing  n  +  t  =  r,  the  second  term  becomes 

oo  r  (-)r  n(k2s2/4)r 
I  I 


(m-1)!  r=0  n=0  (m+n) * ! r!  (r-n)! 


Equating  terms  in  (k2s2/4)r. 


1 

1 

l  (~)n 

(m+r)!  r! 

(m-1)! 

n_g  (m+n)n!  r!  (r-n)! 

1 

_  T  (")n 

m(m+l)(m+2). . 

. . (m+r) 

n_Q  (m+n)r! (r-n)! 

(Eq.  A. 8) 


a  relationship  which  is  surprisingly  difficult  to  establish  by  purely 
algebraic  methods. 


A. 2. 2  Vertical  Array  (y  =  n/2) 

A  similar  simplification  is  introduced  by  setting  y  =  7t/2.  In  this  case, 
in  Eq.  A  4,  the  only  non-zero  term  in  the  summation  over  n  is  that  for 
which  n  =  0,  yielding 


Ii(*  =  n/2) 


I  ^~l<2s2-)  II2m+2i+l(Az) 
i=0  (2i )! 


(Eq.  A. 9) 
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This  is  an  even  simpler  form  than  for  the  horizontal  array.  The  factors 

r  _  .  2 

- - -  -  are  the  terms  in  the  expansion  of  cos(ks). 

(21)! 

The  computation  is  similarly  somewhat  simpler.  As  for  the  horizontal  array 
there  is  a  terminal  condition  that 

(ks)2i+2  (ks)1' 

(21+2)!  ( 2 i  ) ! 

or 

(2i+2)(2i+l)  >  k2s2  , 
which  is  certainly  the  case  if 

(2i+l)  >  ks 
or 

i  >  (ks-l)/2. 


Moreover,  since  terms  alternate  in  sign,  a  sufficient  condition  for  an 
error  not  exceeding  e  is 


(k2s2)i+1 

(2i+2)! 


E2m+2i+l(Az) 

E2m+l(Az) 


<  e 


(Eg.  A. 10) 
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APPENDIX  B 

EVALUATION  OF  THE  INTEGRAL  GIVING  THE  AMBIENT  NOISE  POWER 

IN  THE  OUTPUT  OF  A  CONVENTIONAL  BEAMFORMER 


In  the  main  text  Eq.  4c  includes  the  following  double  integral: 
n/2  2  n 

~  f  ^  cos2m  ■*'P*sinp*exp(-Az*secp)  [f(P)F(Q)]2  d(J>  dp  , 

0  0  (Eq.  B.l) 

where  m  is  a  non- negative  integer, 

Az  is  a  positive  constant, 

frp>  =  sin(Px/2)  F/-n.  sin(Qy/2) 

u  J  P  sin(x/2)  5  w  Q  sin(y/2)  ’ 

with  x  =  ks*sinp*cos<()  , 

y  =  ks(sinp*cosysin<{>  -  sinycosp  +  sinO). 


The  first  stage  is  the  integration  with  regard  to  <|>,  which  involves  only 
the  factor  in  squared  brackets.  To  start  with  we  find  an  expression  for 
f2(P)  in  terms  of  the  weighted  sum  of  a  number  of  terms  of  the  form 
cos(rx).  To  show  that  this  is  possible,  consider  the  trigonometric 
identify: 


sin(Px/2)  -  sin(P-2*x/2)  =  2cos(P-l*x/2)*sin(x/2) 
and  therefore 


sin(P-2)*x/2)  -  sin(P-4*x/2)  =  2cos(P-3*x/2)*sin(x/2) ,  etc. 
terminating  with 

sin(2x/2)  -  sin(O)  =  2cos(x/2)sin(x/2)  [P  even] 
sin(3x/2)  -  sin(x/2)  =  2cos(2x/2)sin(x/2)  [P  odd]. 


Summi ng, 

S1n(Px/2)  =  2[cos(P-T)*x/2)  +  cos(7-3*x/2)  +  ...]  (Eq.  B.2) 

sin  x/2 

the  summation  ending  with  cos(x/2)  [P  even] 

1/2  [P  odd]. 
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On  squaring, 


sin2(Px/2) 

sin2(x/2) 


will  contain  terms  of  the  form 


cos(rx/2)cos(tx/2)  =  %[cos(r+t) -x/2)  +  cos(r-t)-x/2] . 

But  r  and  t  are  either  both  odd  or  both  even,  so  that  ^(r+t)  and  ^(r-t)  are 
both  integers.  It  follows  that 

sin2(P*/2)  x  /r-  r,  -.n 

- *-2 — -  =  I  A  cos(rx)  (Eq.  B.3) 

sin2(x/2)  r=0  r 

[since,  on  squaring  Eq.  B.2,  the  highest  order  term  will  be  cos(P-l)x]. 

To  evaluate  the  coefficients  A^  in  Eq.  B.3,  note  that 

sin2(Px/2)  _  1-cosPx 
sin2(x/2)  1-cosx 


so  that 


P-1 


1-cosPx  =  (1-cosx)  I 
r=0 


A  *cos(rx) 
r  v  J 


P-1 


P-1 


I  A  'cos(rx)  -  A  -cosx  -  h  I  A  [cos(r+lx)  +  cos(r-lx)]. 
r=0  r  0  r=l  r 


Equating  coefficients  of  cos(rx). 


A0  ~  %AX  =  1, 

~  ^2  = 

Ar  -  «r-l  -  %Artl  =  0  2  <  r  <  P-2, 

Ap_i  -  ^Ap_2  =  0, 


It  is  easy  to  show  that  the  solution  of  these  equations  is 


Ar  =  2(P-r)  1  <  r  <  P-1. 
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Simi larly 

s1n2(<fr/2)  ^  B.-cos(ty), 
sin2(y/2)  t=0 

where 

B  =  Q, 
o 

Bt  =  2(Q-t)  1  <  t  <  Q-l  . 

Therefore 

P-1  Q-l 

p2Q2[f(p).f(Q)]  2  i  a  B  cos(rx) *cos(ty) 

r=0  t=0 


=  h  2  2  A  Bt[cos(rx+ty)  +  cos(rx-ty)]. 
r  t  r 


Now  consider  the  integral  of  cos(rx+ty)  with  respect  to  0 . 

rx+ty  =  ks[r-sinP'Cos<t>  +  t(sinp*cosy s i n0  -  sinycosp  +  sin0] 
=  ks[sinp  Vr2+t2cos2y  si n((|)+e)  -  t(siny*cosp  -  sin0)], 


where 


tane  =  r/(t*cosy). 


Therefore 


cos(rx+ty)  =  cos(ks*sinp  A/r2+t2cos2y •  s i n<j)+e ) •  cos [ ks t ( s i ny •  cosp  -  sin0] 


+  sin(ks*sinp  vr2+t2cos2y*sin(})+e)‘sin[kst(sinycosp  -  sin0]. 


But,  in  integrating  over  2 n,  since  the  integrand  is  periodic  with  period 
2 n,  the  phase  term  e  does  not  appear  in  the  result,  and  may  be  therefore 
set  equal  to  zero. 

Further,  in  forming  cos(rx-ty),  it  is  necessary  only  to  change  the  sign  of 
t,  which  will  leave  the  cosine  term  unaltered,  but  will  change  the  sign  of 
the  sine  term.  The  latter  therefore  vanishes  identically  on  adding,  and 

2  n  2  n 

P2q2  f  [f(P)*f(Q)]2d<|>  =  I  I  J  A  Bt*cos[kst(siny*cosp-sin0] 
o  r  t  o  r  x 

cos(A*sinP‘sin<j>)  d4> 
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where  A  =  ks  Vr2+t2cos2y 

=  2n  I  I  A  B  ’cos[kst(sinycosp  -  sine]J  (A-sinp) 
r  t  r  L  0 

=  2n  1  I  A  B.J  (A*si np) [cos(kst*si n0)cos(kst*si nycosp) 
r  t  r  1  ° 

+  sin(kst«sin0)sin(kst-siny cosp)] , 


and  therefore 

P2Q2  P-1  Q-l  ti/2 

-  1  =  2  2  /  A  B. [cosp-cos(qcosp)+sinp*sin(q*cosp)]J  (A*sinp) 

2/t  2  r=0  t=0  o  r  r  o 

x  cos2m  ^‘sinp-expC-Az'secp)  dp  , 

where 

p  =  kst*sin0 
q  =  kst-siny 
A  =  ks  7  r2+t2cos2y  . 


Now  the  integral  involving  cos(q-cosp)  is  of  the  same  form  as  that  for 
evaluating  the  integral  Ix  discussed  in  Appendix  A,  and  we  may  write 
immediately  the  first  term  as  (see  Eq.  A. 4  of  App.  A). 


oo  po 
2  2 
i=0  n=0 


P-1  Q-l 

I  I  A  B. *cosp 
r=0  t=0  r 


(-d2)1 

(21)! 


Jn(M  E2m+2i+2n+l 


(Az). 


Moreover,  in  forming  the  second  integral,  it  is  sufficient  to  note  that 
s  i  nee 


cos(q*cosp)  = 


y  (-q2cos2p)1 

1=0  <21>! 


and 


sin(q*cosp) 


00 

I 

i=0 


?i+1 

(q-cosp) 

}  (2i+l)! 


=  q*cosp 


t  (~q2cos2P)1 

1=0  <2i+1>! 


this  term  derives  immediately  from  the  first  by  writing  2m+l  for  2m  and 
writing  (2i+l)  for  (2i). 
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Thus  finally, 

,2  00  00  P-1  Q-l 


P2QZ 

2  n 


*>  =io  ^  J"(X)  '  T^E2^n+2in<Az> 


where 


A  =  P 
o 

Ar  =  2(P-r)  [r>0] 


k  =  ks  V  r2+t2cos2y 
p  =  kst*sin0 
q  =  kst-siny 


si np*q 


2i+l 


(2i+l)! 


E2m+2n+2i+2<-AZ')  ’ 


(Eq.  B. 4) 


B  =  Q 
o  x 

Bt  =  2(Q-t)  [t>0] 


This  formidable  quadruple  integral  is  the  formal  solution  of  the  double 
integral.  The  general  form  (Eq.  B.4)  has  been  derived  for  the  sake  of 
completeness  although  it  does  not  appear  to  be  inviting  for  numerical 
computation.  It  is  not,  however,  impossible.  Comparing  with  the 
equivalent  form  for  Ix  (Eq.  A. 4  of  App.  A),  it  is  evident  that,  for  given  n 
and  i,  Eq.  A. 4  requires  only  one  Bessel  function  and  one  exponential 
integral  to  be  computed.  In  the  present  integral  it  is  immediately  evident 
that  PQ  Bessel  functions  will  need  to  be  evaluated,  two  exponential 

integrals,  and  2PQ  coefficients  of  the  form  A  B.  q21  or  A  B  q2  +1.  For 
PQ  large  this  would  be  a  formidable  undertaking.  r 

However,  it  is  probably  much  more  important  to  look  at  the  special  cases  of 
y  =  0  and  y  =  n/2. 

Making  y  =  0  does  yield  notable  simplification.  Since  q  =  0,  the  only 
surviving  term  in  the  summation  over  i  is  the  first  term  in  Eq.  B.4  for 
i  =0,  so  that 


oo  P-1  Q-l 

I2  (y=0)  “ - 2  2  2  A  B  *cosp 

P2Q2  n=0  r=0  t=0 

where 

Ax  =  ks  J  r2+t2  p  =  kst-sine  . 


ax/2)' 


Jn(Al)E2m+2n+l(Az): 


(Eq.  B. 5) 


Comparing  with  the  equivalent  form  for  Ix  (Eq.  A. 5  of  App.  A),  Eq.  B.5 
consists  of  a  weighted  sum  of  PQ  simpler  forms,  so  that  a  computer 
programme  designed  for  Ix  could  be  used  as  a  subroutine  for  I2.  This  would 


37 


SACLANTCEN  SM-159 


be  slightly  wasteful  in  the  computation  of  the  Bessel  Functions  since, 
because  of  the  symmetry  of  in  r  and  t,  many  may  turn  up  twice.  It  is 
easy  to  show  in  fact,  if  R  is  the  lesser  of  P  and  Q,  the  number  of 
different  Bessel  functions  to  be  calculated  is  PQ-Js( R-l) .  For  P  =  Q,  for 
example,  this  becomes  SgPCP+l),  a  saving  of  nearly  50%  for  P  large. 
However,  the  complication  in  programming  for  these  replicas  may  offset  such 
a  savi ng. 

Making  y  =  n/2  does  not  simplify  quite  so  neatly.  We  now  retain  the 
quadruple  summation,  so  that  I2 ( Y  =  n/2)  has  the  form  of  Eq.  B.4,  in  which 

y2=  ksr  q  =  kst. 


However,  \2 

(W  2)n  j 

n 


does  not  involve  t,  so  there  are  now  only  P  terms  of  the  form 
(\2)  to  be  computed. 
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